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Abstract. This paper considers the valuation of exotic path-dependent 
options in Levy models, in particular options on the supremum and the 
infimum of the asset price process. Using the Wiener-Hopf factoriza- 
tion, we derive expressions for the analytically extended characteristic 
function of the supremum and the infimum of a Levy process. Combined 
with general results on Fourier methods for option pricing, we provide 
formulas for the valuation of one-touch options, lookback options and 
equity default swaps in Levy models. 



1. Introduction 

The ever-increasing sophistication of derivative products offered by finan- 
cial institutions, together with the failure of traditional Gaussian models to 
describe the dynamics in the markets, has lead to a quest for more realistic 
and flexible models. In fact one of the lessons from the current financial crisis 
is the following: the Gaussian copula model is inappropriate to describe the 
interdependence between the tails of asset returns because, among other pit- 
falls, the tail dependence coefficient is always zero; hence, this model cannot 
capture systemic risk. 

In the search for appropriate alternatives, Levy processes are playing a 
leading role, either as models for financial assets themselves, or as building 
blocks for models, e.g. in Levy-driven stochastic volatility models or in affine 
models. The field of Levy processes has become popular in modern mathe- 
matical finance, and the interest from academics and practitioners has led 
to inspiring and challenging questions. 

Levy processes are attractive for applications in mathematical finance be- 
cause they can describe some of the observed phenomena in the markets in a 
rather adequate way. This is due to the fact that their sample paths may have 
jumps and the generated distributions can be heavy-tailed and skewed. An- 
other important improvement concerns the famous smile effect. See Eberlein 



and Keller (119!)") i and Eberlein and Prause (2002) for an extensive empirical 



justification of the non-Gaussianity of asset returns and the appropriateness 
of (generalized hyperbolic) Levy processes. For an overview of the appli- 
cation of Levy processes in finance the interested reader is referred to the 
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textbooks of Cont and Tankov (2004) , Schoutens (2003) as well as the col- 



lection edited by Kyprianou et al. (2005) There are, of course, several text 



books dealing with the theory of Levy processes; we mention Bertoin (1996) 



Sato (|1999p, Applebaum (2004) and Kyprianou (2006) while the collection 



by Barndorff-Nielsen et al. (2001) contains an overview of the application of 



Levy processes in different areas of research, such as quantum field theory 
and turbulence. 

The application of Levy processes in financial modeling, in particular for 
the pricing and hedging of derivatives, has led to new challenges of both an- 
alytical and numerical nature. In Levy models simple closed form valuation 
formulas are typically not available even for plain vanilla European options, 
let alone for exotic path-dependent options. The numerical methods which 
have been developed in the classical Gaussian framework lead to completely 
new challenges in the context of Levy driven models. These numerical meth- 
ods can be classified roughly in three areas: probabilistic numerical methods 
(Monte Carlo methods), deterministic numerical methods (PIDE methods), 
and Fourier transform methods; for an excellent survey of these methods, 



their applicability and limitations, we refer to Hilber et al. (2009) 



This paper focuses on the application of Fourier transform methods for 
the valuation of exotic path-dependent options, in particular options de- 
pending on the supremum and the infimum of Levy processes. The bulk 
of the literature on this latter topic focuses on the numerical aspects. Our 
focus is on the analytical aspects. More specifically, we show first that the 
Wiener-Hopf factorization of a Levy process possesses an analytic exten- 
sion, and then we prove that the Wiener-Hopf factorization (viewed as a 
Laplace transform in time) can be inverted. These results allow us to derive 
expressions for the extended characteristic function of the supremum and 
the infimum of a Levy process. This latter result, combined with general re- 
sults on option pricing by Fourier methods (cf. lEberlein et al. 2010|) , allows 
us to derive pricing formulas for lookback options, one-touch options and 
equity default swaps in Levy models. 

Let us briefly comment on some papers where the Wiener-Hopf factoriza- 
tion is used to price exotic options in Levy models. Boyarchenko and Leven- 
dorskh (2002a]) derive valuation formulas for barrier and one-touch options 
for driving Levy processes that belong to the class of so-called "regular Levy 
processes of exponential type" (RLPE); cf. also the book by Boyarchenko 
and Levendorskh (!2 002bp . The results of these authors are based on the the- 
ory of pseudodifferential operators. The numerics of this approach is pushed 



further in Kudryavtsev and Levendorskh (2006, 2009). Avram et al. (2004) 



Asmussen et al. (2004) , Kyprianou and Pistorius (2003) , Alili and Kypri 



anou ()2005p . and Levendorskh et al. (2005) consider the valuation of Amer- 
ican and Russian options, either on a finite or an infinite time horizon. 
Jeannin and Pistorius (|2009j) develop methods for the computation of prices 
and Greeks for various Levy models. Central in their argumentation is the 
approximation of different Levy models by the class of "generalized hyper- 
exponential Levy models", which have a tractable Wiener-Hopf factoriza- 



tion. The same approach is also applied in Asmussen et al. (2007) for the 
pricing of equity default swaps in Levy models. 
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The major open challenge in this field is the development of analytical 
expressions for the Wiener-Hopf factors for general Levy processes. In a 



remarkable recent development, Hubalek and Kyprianou (2010) generate a 



family of spectrally negative Levy processes with tractable Wiener-Hopf 
factors, using results from potential theory for subordinators. These re- 



sults were later extended in Kyprianou and Rivero (2008) and applied to 



problems in actuarial mathematics in Kyprianou et al. (2009) Moreover, in 



two very recent papers IKuznetsovl (|2010|. I2009P introduces special families 
of Levy processes such that the Wiener-Hopf factors can be computed as 
infinite products over the roots of certain transcendental equations. These 
families include processes with behavior similar to the CGMY process, while 
the author shows that the numerical computation of the infinite products 
can be performed quite efficiently. 

This paper is structured as follows: in section [21 we briefly review Levy 
processes and prove the analyticity of the characteristic function of the 
supremum. In section [3l we review the Wiener-Hopf factorization, prove 
its analytic extension and invert it in time. In section U we present some 
examples of popular Levy models and comment on the continuity of their 
laws. Finally, in section [SJ we derive valuation formulas for lookback and 
one-touch options as well as for equity default swaps. 

Important Remark. This paper is intimately tied to, and intended to be 



read together with, the companion paper |Eberlein, Glau, and Papapantoleon 
(|2010p . which will be abbreviated EGP in the sequel. In particular, we will 
make heavy use of the notation and results from that paper. 



2. Levy processes 

We start by fixing the notation that will be used throughout the paper 
and providing some estimates on the exponential moments of a Levy process. 
Then, we prove the analytic extension of the characteristic function of the 
supremum and the infimum of a Levy process, sampled either at a fixed time 
or at an independent, exponentially distributed time. 



2.1. Notation. Let ¥> = (O, J 7 , F,P) be a complete stochastic basis in the 
sense of Jacod and Shiryaev (2003, 1.1.3) where T = Ft, < T < oo 
and F = (J r t)o<t<T- Let L = (L t )o<t<T be a Levy process on this stochastic 
basis, i.e. L is a semimartingale with independent and stationary increments 
(PUS), and Lq = a.s. We denote the triplet of predictable characteristics 
of L by (B,C,v) and the triplet of local characteristics by (6, c, A); using 



Jacod and Shiryaev (2003, II. 4. 20) the two triplets are related via 



B t {oj) = bt, Ct(uj) = ct, v(u};dt,dx) = A(dx)di. 
We assume that the following condition is in force. 
Assumption (EM). There exists a constant M > 1 such that 



e ux X(dx) < oo, 



Vu e [-M,M]. 



{\x\>l} 
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The triplet of predictable characteristics of a PUS determines the law of 
the random variables; more specifically, for a Levy process we know from 
the Levy-Khintchine formula that 

E[e iuLt ] = exp (t • K(iu)) , (2.1) 

for all t G [0, T] and all where the cumulant generating function is 

u 2 f 

k(u) = ub+—c+ / {e ux - 1 - ux)X(dx). (2.2) 

R 

Assumption (EM) entails that the Levy process L is a special and expo- 
nentially special semimartingale, hence the use of a truncation function can 



be and has been omitted. Applying Theorem 25.3 in Sato (1999) we get that 



E [e uLt ] < oo, Vu G [— M, M], Vt G [0,T]. 

Recall that for any stochastic process X we denote by X the supremum 
and by X_ the infimum process of X respectively. 

In the sequel, we will provide the proofs of the results for the supremum 
process. The proofs for the infimum process can be derived analogously or 
using the duality between the supremum and the infimum process; see the 
following remark. 

Remark 2.1. Let L be a Levy process with local characteristics (b, c, A). 
The dual of the Levy process L defined by L' := — L, has the triplet of local 
characteristics (6', c', A') where b' = —b, c' = c and 1a{%) * A' = 1a{—x) * A, 
A G B(R\{0}). Moreover, we have that 

L t = inf L s = - sup (-L s ) = -L' t . 

0<s<t 0<s<t 

2.2. Analytic extension, fixed time case. In this section, we establish 
the existence of an analytic extension of the characteristic function of the 
supremum and the infimum of a Levy process, and derive explicit bounds 
for the exponential moments of the supremum and infimum process. 

The next lemma endows us with a link between the existence of exponen- 
tial moments of a measure g and the analytic extension of the characteristic 
function g. 

Lemma 2.2. Let g be a measure on the space (M,/3(1R)). If J e ux g(dx) < oo 
for all u £ [—a, b] with a,b > 0, then the characteristic function g has an 
extension that is continuous on (—00,00) x i[— b, a] and is analytic in the 
interior of the strip, (—00,00) x i{—b,a). Moreover g(u) = J e lux g(dx) for 
all u G C with Q(u) G [—6, a]. 

Proof. The function u 1— >■ e lux clearly extends to an entire function and the 
extension 

g{u) := J e lux g(dx) (116C with 3?(tt) G [-6, a]) 
is well-defined since 

\J.ux\ — Q(u)x ^ „-aii I Jyx-t 

|e |-e <e l{x<o} + e i -{x>o} =■ h{x), 
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for u G C with 9(n) G [—6, a], and we have that /i G L 1 (g) by assump- 
tion. Moreover, Lebesgue's dominated convergence theorem yields that this 
extension is continuous. 

We will prove the analyticity of gin (—00, 00) x i(—b, a) using the theorem 
of Morera (cf. for example Theorem 10.17 in lRudin 1987ft . Let 7 be a triangle 
in the open set (—00,00) x i(—b,a); the theorems of Fubini and Cauchy 
immediately yield 

d-y d-y d-y 

as « 4 e lux is analytic for every fixed x G R. Then, the analyticity of 
g follows from Morera's theorem. For a justification of the application of 
Fubini's theorem it is enough to note that 



du g(dx) < J J h(x)du g{dx) = £(7) J h(x)g(dx) < 00, 

d-y d-y 

where £(7) denotes the length of the curve dj. □ 

Lemma 2.3. Let Y be a Levy process and a special semimartingale with 
E[Y t ] = for some, and hence for every, t > 0. Then 

E[e Y *] <8£[el y *l], 

where Y* = sup < s < t \Y S \. 

Proof. Using that - ' / is positive for every n > and the monotone con- 
vergence theorem, we get 



n\ 1 — ' n\ 

n=0 n=0 



Now, Remark 25.19 in |Sato (1999)] yields 

E(Y t *) n < 8E\Y t \ n , for every n > 1, 
while for n = the inequality holds trivially. Hence, we get 

(V*\n 00 |v|n 00 IV" 



y E ifi^< 8 y E ^L = 8jB y JUL ^u^h □ 

^— ' n! ^— ' n! ^— ' n! L J 

n=0 n=0 n=0 

Next, notice that under assumption (EM) we have that 

J \e Mx - 1 - Mx\ \{dx) < 00 and y |e" Mx - 1 + Mx\ \{dx) < 00. 

R R 
Let us introduce the following notation: 

a(M) := M\b\ + ^cM 2 + J \e Mx - 1 - Mx\ X(dx) (2.3) 



and 



a(M) := M\b\ + ^cM 2 + J \e~ Mx - 1 + Mx\ X(dx). (2.4) 
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Lemma 2.4. Let L = (Lt)o<t<T be a Levy process that satisfies assumption 
(EM). Then we have the following estimates 

E[e uTt ] < E[e MTt ] < 8Q(t, M) < oo (u < M), 

and 

E[e~ u ^] < E[e~ M ^] < 8Q(t, M) < oo (u < M), 
where Q(t, M) := e m ^ + e ^ M \ 
Proof. For u<Mwe have 

e uIl <e Mli , 

since L t = sup 0<s<t L s is nonnegative. Further notice that 

~L t = sup [bs + ■s/cW s + Lf\ < sup [s/cW s + L d s \+ sup [6s], 

0<s<t 0<s<t 0<s<t 

where Lt = bt + \fcWt + L d denotes the canonical decomposition of L, with 
Brownian motion W and a purely discontinuous martingale L d = x*(fi — u). 
Let us further denote by 

Y s := JdW 8 + L d s . 

The process Y is not only a martingale but also a Levy process and a special 
semimartingale with local characteristics (0, c, A). We have 

L t < sup Y a + \b\t < Y* + \b\t, 

0<s<t 

hence we get that 

E[e MTt ] < E [e M ^+^] = e M W E[e MY *] < 8^^E[e M ^] , (2.5) 

using Lemma 12.31 for the special semimartingale Z := MY, which is a Levy 
process satisfying E[Zt] = for every < t < T. 
Now it is sufficient to notice that 

E[ e M \ Y ^} < E[e MYt ] + E[e~ MYt ] , (2.6) 



where Theorem 25.17 in Sato (1999) yields 



E[e MYt ] = exp (t^- + t J {e Mx -l-Mx) A(dx)) 

R 

< e (HM)-M\b\)t. (2 _ 7) 

similarly, 

E[e~ MY ] < e fe(M)-M| b |)t_ (2.8) 
Summarizing, we can conclude from (|2.5p - ([2.8p that 

E[e MTt ] < 8e A/|fe|i (e (s(M) " M|6|)< + e (s(M) ~ M|6|)t ) 

as well as 

E[e- M ±<] < 8(e^ M )* + ^1 . □ 
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A corollary of these results is the existence of an analytic continuation for 
the characteristic function ipj^ of the supremum, resp. (pL t of the infimum, 
of a Levy process. 

Corollary 2.5. Let L be a Levy process that satisfies assumption (EM). 
Then, the characteristic function (pj^ of Lf, resp. ifL t of L t , possesses a 
continuous extension 

= J e izx P Tt (dx), resp. ^(z) = J P kt (dx) , 

R R 
to the half-plane z £ {z £ C : —M < Qz}, resp. z £ {z £ C : < M}, 
that is analytic in the interior of the half-plane {z € C : — M < $sz}, resp. 
{z£C:Qz< M}. 

Proof. This is a direct consequence of Lemmata 12.21 and 12.41 □ 

Remark 2.6. One could derive the statement of Corollary 12.51 using the 
submultiplicativity of the exponential function and Theorem 25.18 in Sato 



(|1999p . see Lemma 5 in Kyprianou and Surya (2005) However, we will need 



the estimates of Lemma 12.41 in the following sections. 

2.3. Analytic extension, exponential time case. The next step is to es- 
tablish a relationship between the (analytic extension of the) characteristic 
function of the supremum, resp. infimum, at a fixed time and at an inde- 
pendent and exponentially distributed time. Independent exponential times 
play a fundamental role in the fluctuation theory of Levy processes, since 
they enjoy a property similar to infinity: the time left after an exponential 
time is again exponentially distributed. 

Let 9 denote an exponentially distributed random variable with parameter 
q > 0, independent of the Levy process L. We denote by Lg, resp. Lg, the 
supremum, resp. infimum, process of L sampled at 9, that is 

Lg = sup L u and Lg = inf L u . 

0<u<9 0<u<6 

Lemma 2.7. Let L = (Lt)o<t<T be o> Levy process that satisfies assumption 
(EM), and let 9 ~ Exp(q) be independent of the process L. 

If q > a(M) Va(M), then the characteristic function (pj^ of Lg possesses 
a continuous extension 

oo 

= I * IZX P-lM x ) = 1 1 e-^e^d* (2.9) 

R 

to the half-plane z £ {z E C : — M < Qz}, that is analytic in the interior of 
the half-plane {z £ C : — M < Qz}. 

Lf q > a(M) Va(M), then the characteristic function (fL of L_g possesses 
a continuous extension 

oo 

VhM = j * izx PlM x ) = 1 J e~ qt E[e lz ^]dt (2.10) 

R 
to the half-plane z £ {z £ C : Qz < M}, that is analytic in the interior of 
the half-plane {z £ C : Qz < M}. 
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Proof. We have that 

oo oo oo 

E[e uTe ) = J J e ux qe- gt P Tt (dx)dt = J E [e^] qe^dt, 



and, for q > a(M) V a(M), by Lemma l2~4l we get 



;[ 9 -a(M)] 



dt + q J e * 
o 



[ 9 -a(M)] 



di < oo; 



E[e MLt ]qe- qt dt <8iq J e" 
^ 

hence, for u < M, we have 

£[e nX "] < £[e Ml «] < oo (g > 5(M)V 2 (M)). (2.11) 

Inequality (|2.1ip . together with Lemma |2.2|, implies that the characteristic 
function <p-^ has a continuous extension to the half-plane {z G C : — M < 
that is analytic in {z G C : — M < and is given by 



E[ 



jzLg] 



for every z G C with 9z > — M. Furthermore Fubini's theorem yields 

oo oo oo 

£[ e ^ fl ] = /" /" e i2 ^ e - ,? *P i;t (dx)dt = g f e- qt E[e izTt ]dt. 





The application of Fubini's theorem is justified since, for Qz > — M and 

q > a(M) Va(M), we have 



Ml 



-AzLa 



E[e~^ z)Le ] < E[e ML °] < oo 



by inequality (|2.1ip . Similarly, we prove the assertion for the infimum. □ 

3. The Wiener-Hopf factorization 

We first provide a statement and brief description of the Wiener-Hopf 
factorization of a Levy process, and then show that the Wiener-Hopf fac- 
torization holds true for the analytically extended characteristic functions. 
Next, we invert the Wiener-Hopf factorization, and derive an expression for 
the (analytically extended) characteristic function of the supremum, resp. 
infimum, of a Levy process in terms of the Wiener-Hopf factors. 

3.1. Analyticity. Fluctuation identities for Levy processes originate from 
analogous results for random walks, first derived using combinatorial meth- 
ods, see e.g. Spitzer (1964) or Feller (1971) Bingham (1975) used this 



discrete-time skeleton to prove results for Levy processes; the same approach 



is followed in the book of |Sato (1999)} IGreenwood and Pitmanl (|198Uall98Ubp 
proved these results for random walks and Levy processes using excursion 



theory; see also the books of Bertoin (1996) and Kyprianou (2006) 

The Wiener-Hopf factorizatior^ serves as a common reference to a mul- 
titude of statements in the fluctuation theory for Levy processes, regarding 
the distributional decomposition of the excursions of a Levy process sam- 
pled at an independent and exponentially distributed time. The following 



^The historical reasons leading to the adoption of the terminology "Wiener-Hopf" are 



outlined in section 6.6 in Kyprianou (2006) 
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statement relates the characteristic function of the supremum, the infimum, 
and the Levy process itself. Let L be a Levy process and 9 an independent, 
exponentially distributed time with parameter q; then we have that 

E[e izLe ] = E[e izTe ]E[e izL e] 

or equivalently, 



here k denotes the cumulant generating function of L±, cf. (|2.2p . and ip~£~, 
ip~ denote the so-called Wiener-Hopf factors. 

In the sequel, we will make use of the Wiener-Hopf factorization as stated 



in the beautiful book of Kyprianou (2006) , and prove the analytic extension 



of the Wiener-Hopf factors to the open half-plane {zdC: Qz > — M}. 
Recall the definitions of (|2.3p and (|2.4p . and let us denote by 



a*(M) := max|a(M),a(M)|. 



Theorem 3.1 (Wiener-Hopf factorization). Let L be a Levy process that 
satisfies assumption (EM) (and is not a compound Poisson process). The 
Laplace transform of Lg, resp. Lq, at an independent and exponentially dis- 
tributed time 9,9^ Exp(g), with q > a*(M), can be identified from the 
Wiener-Hopf factorization of L via 

oo 

E[e-M]=[ q E[e-K]e-«K=^ (3.1) 



o 



and 



oo 

£[ ^ ]= | gB[e ^ ]e -^ = || (3.2) 
o 

for (3 £ {/3 £ C : 9i(/3) > — M}. The Laplace exponent of the ascending, 
resp. descending, ladder process 7t(a,/3), resp. tk{a,f3), for a > a*(M) and 
k,k> 0, has an analytic extension to (3 £ {/3 £ C : 9?(/3) > — M} and is 
given by 

«(a,/9) = &exp [ / / (e"* -e" Q *- /3a; )ip Lt (dx)di ), (3.3) 




and 



K(a,/3) = fcexp ^ ^ (e~* - e~ at+l3x )^P Lt (dx)dt 



(3.4) 



(-oo,0) 



Remark 3.2. Note that the Wiener-Hopf factors and (p~ are related 
to the Laplace exponents of the ascending and descending ladder process k 
and k via 

^ (z/3) = ^y and v < { - %0) = tifc0)' (3 - 5) 
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We will prepare the proof of this theorem with an intermediate lemma. 
Let us denote the positive part by a + := max{a, 0}. 

Lemma 3.3. Let L be a Levy process that satisfies assumption (EM). For 
q > k(M) + the maps 

oo 

z^J j {l-e lzx )P Lt (dx)^dt (3.6) 

(0,oo) 

and 

oo 

z^ j j (e^ -e- gt+izx )P Lt (dx)^dt (3.7) 

(0,oo) 

are well defined and analytic in the open half plane {zdC: ^( z ) > — M\. 

Proof. We will show that for every compact subset K C {z € C : Q(z) > 
— M}, there is a constant C = C(K) > such that 

oo 

J j \e izx -l\P Lt (dx)^-dt<C(K), (3.8) 

(0,oo) 

for every z £ K. Then, applying Lebesgue's dominated convergence theorem 



yields the continuity of the function 

oo ^ 

z^ J J (e izx -l)P Lt (dx)^dt 

(0,oo) 

inside the half-plane {z S C : $s(z) > — M}. Moreover, let 7 be an arbitrary 
triangle inside {z G C : $s(z) > — M}; the theorems of Fubini and Cauchy 
yield 



J (e izx -l)P Lt (dx)^-dtdz 

d-y (0,oo) 

00 

= // J (e lzx -l)dzP Lt (dx)^dt = 0. (3.9) 

(0,00) 97 

Hence, applying Morera's theorem yields the analyticity of ()3.6p in the open 
half-plane {z £ C : $S(z) > —M}. 

The assertion for the second map immediately follows from the identity 

(e~* - e-*****)*- 1 = (1 - e^e -9 ** -1 + (e - * - e -9 *)* - 1 

and the integr ability of the second part, since 

00 

|e~* - e _9f |t _1 dt < 00 
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and 



j \e~' - e""|( _1 di = J \e' { «~ 1] - l|e""( _1 dt < C\g - 1\ J e~"dt < 



OO, 



with C > 1, for e > small enough. 

To show estimation (|3.8p . we choose a constant k = k(K) > only 
depending on the compact set K, such that \z\ < k for every z G K, and we 
write 



(0,oo) 



J \e izx - l\P Lt (dx) 

J \e izx - l\P Lt (dx) + y |e^-l|P it (dx) 

(o,i/fc] 

< J \zx\P Lt (dx)+ j \e lzx \P Lt (dx)+ j P Lt (dx). (3.10) 



(o,i/fc] 



(l/fc,oo) 



(l/fc,oo) 



Using inequality (30.13) of Lemma 30.3 in Sato (1999) we can deduce 
j \zx\P Lt (dx) < k I \x\P Lt (dx) < kE[\L t \l {{Ltl < 1/k} ] < diK)^ 2 , 



(o,i/fe] 



(0,1/fc] 



with a constant C\{K) that depends only on the compact set K. Similarly, 
using inequality (30.10) in Sato (1999) we can estimate the last term of 

j P Lt (dx) = P{{L t > 1/k}) < P({\Lt\ > 1/fc}) < C 2 (K)t, 

(l/fe,oo) 

with a constant C2(K) that depends only on the compact set K. In order to 
estimate the second term of inequality (|3.10p , let us note that we may choose 
e > small enough, such that for every z G K, we have — < M' < M 
with M' := M(l - e), and we get 

J \e iz *\P Lt (dx)<E[e M ' L n {lLtl>1/k} ). 

(l/fe,oo) 

Applying Holder's inequality withp := j-^ and q := -, together with Lemma 
30.3 in|Sato (1999)1 yields 



i/p 



E[e M ' Lt l {lLtl>1/k} ] < (E[eO M ' L *]y(P({\L t \ > 1/k} 
< C 3 (K)t e e {1 - e)K(M)t . 

Altogether we have 

J \e lzx - l\P Lt {dx) < Ci(i^)t 1/2 + C 2 (K)t + C 3 (K)t e e^ K 



(AO* 



(0,oo) 
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with positive constants C\(K), C^-fQ and C%{K) that only depend on the 
compact set K. As q > (1 — e)(n(M)) + , we can conclude (|3.8p . which com- 
pletes the proof. □ 

Proof of Theorem For (3 G C with > the assertion follows directly 



from Theorem 6.16 (ii) and (iii) in Kyprianou (2006) 



From Lemma 12.71 we know that for q > a*(M) the function 

/3*<p Ia (iP)=E[e-f SL °] 
has an analytic extension to the half-plane 

{P G C : »03) > — M}, 
whereas Lemma 13.31 yields that if q > a*(M), the mapping 

has an analytic extension to the half-plane 

{P G C : »(/3) > —M}, 

while identity (|3.3|) still holds for this extension. The identity theorem for 
holomorphic functions yields that equation (|3.ip holds for every {/3 G C : 
&(/3) > — M} if q > a*(M). The proof for equations fl£2H and (J33D follows 
along the same lines. □ 

Remark 3.4. Note that, by analogous arguments, we can prove that the 
Laplace exponent of the ascending, resp. descending, ladder process ~K(a,P), 
resp. k(cx,P), has an analytic extension to a G {a G C : 5R(a) > a*(M)}, 
which is given by (|3.3p . resp. (|3.4p . 



3.2. Inversion. The next step is to invert the Laplace transform in the 
Wiener-Hopf factorization in order to recover the characteristic function 
of Lt, at a fixed time t. Let us mention that although the Wiener-Hopf 
factorization and the characteristic function of Lg are discussed in several 
textbooks, the extended characteristic function of Lt at a fixed time has not 
been studied in the literature before. 

The main result is Theorem 13. 6\ which will make use of the following 
auxiliary lemma. 

Lemma 3.5. The maps 1 1— > ^'[e - ' 3 ^'] and 1 1— > Ele?— *] are continuous for 
all p G C with ftp G [-M,oo). 

Proof. Since the Levy process L is right continuous, stochastically continu- 
ous and L is an increasing process, we get that L s /• Lt a.s. as s — > t. 
As L s > we have 

and we may apply the dominated convergence theorem to get 

S[e _/3ls ] -> J E[ e - /3Xt ] as s -> t, 

for every /3 G C with 3?(/3) > — M. Analogously, taking into account that 
|e^— 3 1 < e~ M — B for > — M, the dominated convergence theorem yields 
the continuity of the second map. □ 
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Theorem 3.6. Let L be a Levy process that satisfies assumption (EM) (and 
is not a compound Poisson process). The Laplace transform of L t and L t at 
a fixed time t, t £ [0, T], is given by 



1 7 e '( y+ ™) k(Y + iv, 0) , 
lim — / — : — : — —dv, 



and 



A-^oo 2ir J Y + iv k(Y + iv, j3) 
-A 



I f e t(Y+iv) ,y + i v 0) , 
lim — / — ~ dv, 



(3.11) 



A^oo 2ir y + iv k(Y + iv, -/3) 



(3.12) 



for (3 eC with K/3 G (-M,oo) and Y,Y > a*(M). 

Proof. Theorem 13 .1\ together with equation (|3.ip . immediately yield 

oo 

1 K(q, 0) 



(3.13) 



for p G C with »(/9) > —M and q > a*(M). 

In order to deduce that we can invert this Laplace transform, we want 



to verify the assumptions of Satz 4.4.3 in Doetsch (1950) for the real and 
imaginary part of 1 1— > E[e~@ Lt ~\ . From the proof of Lemma 12.71 we get that 

oo oo 

J e-^\E[e^] \dt < J e- qt E[e-^ Tt ]dt < oo; 


this yields the required integrability, i.e. absolute convergence, of 



-gt 



Q(E[e^ Lt ])\dt and 



-qt 



&{E[eP Lt ])\dt, 



for q > a*(M). Further the real and imaginary part of 1 1— > E[e~P Lt ~\ are of 
bounded variation for f5 G C with 3i/3 G (— M, oo). 

Let us verify this assertion for the imaginary part, for — M < < 

and < 0. We have that 



-PL, 



iE[ 



sm 



We can decompose sin(x) = f{x) — g{x), where / and g are increasing 
functions with /(0) = g(0) = 0, and |/(x)| < x and \g(x)\ < x. It follows 
that 

sin ( - 9(/3)Z t )e-«< = /( - 9(/3)I t )e- R ^ - g{ - ^)L t )e~W\ 
where both terms are increasing in time and are integrable, since 



E 



-9t(P)L t 



< p{f3)\E 



\Lt\e 



< const -E[e MLt ] < 



oo, 



for h = g and h = f. The assertion for the other parts follows similarly. 
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Now, using the continuity of the map t i— >■ £7[e ^ Lt ], cf. Lemma 13.51 we 



may apply Satz 4.4.3 in Doetsch (1950) to invert this Laplace transform; 
that is, to conclude that 

Y+ioo 
Y— ioo 

= lim — / — )- -K^v, 3.14 

A^oo 2vr ./ y + iu 7e(Y + it;, /3) ' v ; 

-A 

for all (3 G C with G (-M, oo) and for every y > a*{M). The proof for 
the infimum follows along the same lines. □ 

4. Levy processes: examples and properties 

We first state some conditions for the continuity of the law of a Levy 
process, and the continuity of the law of the supremum of a Levy process. 
Then, we describe the most popular Levy models for financial applications, 
and comment on their path and moment properties which are relevant for 
the application of Fourier transform valuation formulas. 

4.1. Continuity properties. The valuation theorem for discontinuous pay- 
off functions (Theorem 2.7 in EGP), and the analysis of the properties of 
discontinuous payoff functions (Examples 5.2, 5.3 and 5.4 in EGP), show 
that if the measure of the underlying random variable does not have atoms, 
then the valuation formula is valid as a pointwise limit. Thus, we present 
sufficient conditions for the continuity of the law of a Levy process and its 
supremum, and discuss these conditions for certain popular examples. 

Statement 4.1. Let L be a Levy process with triplet (b, c, A). Then, The- 



orem 27.4 in Sato (1999) yields that the law Pi t , t G [0, T], is atomless iff L 
is a process of infinite variation or infinite activity. In other words, if one of 
the following conditions holds true: 

(a) : c / or ! { \ x \< 1} |x|A(dx) = oo; 

(b) : c = 0, A(R) = oo and L ^j, |x|A(dx) < oo. 

Statement 4.2. Let L be a Levy process and assume that 

(a) : L has infinite variation, or 

(b) : L has infinite activity and is regular upwards. Regular upwards 
means that P(tq = 0) = 1 where tq := inf{i > : Lt{uj) > 0}. 



Then, Lemma 49.3 in Sato (1999) yields that Lt has a continuous distribu- 
tion for every t G [0, T]. The statement for the infimum of a Levy process is 
analogous. 

4.2. Examples. Next, we describe the most popular Levy processes for ap- 
plications in mathematical finance, namely the generalized hyperbolic (GH) 
process, the CGMY process and the Meixner process. We present their char- 
acteristic functions, which are essential for the application of Fourier trans- 
form methods for option pricing, and its domain of definition. We also discuss 
their path properties which are relevant for option pricing. For an interesting 
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survey on the path properties of Levy processes we refer to Kyprianou and 
Loeffen ([2005b . 

Example 4.3 (GH model). Let H = (H t )o<t<T be a generalized hyper- 
bolic process with C(Hi) = GH(A, a, (3, 8, (J), cf. |Eberlein (2001, p. 321)1 or 



Eberlein and Prause (2002) The characteristic function of Hi is 



a 



2 



P 2 \ ^ K\ (S^/a 2 - (/3 + in) 2 ) 



where K\ denotes the Bessel function of the third kind with index A (cf. 
Abramowitz and Stegun 1968[ ); the moment generating function exists for 



u G (—a — f3,a — (3). The sample paths of a generalized hyperbolic Levy 
process have infinite variation. Thus, by Statements 14.11 and 14.21 we can 
deduce that the laws of both a GH Levy process and its supremum do not 
have atoms. 

The class of generalized hyperbolic distributions is not closed under con- 
volution, hence the distribution of Ht is no longer a generalized hyperbolic 
one. Nevertheless, the characteristic function of C{Ht) is given explicitly by 

<PH t (v) = (<PH 1 (u)) t ■ 

A class closed under certain convolutions is the class of normal inverse 
Gaussian distributions, where A = — ^; cf. Barndorff-Nielsen (1998) . In that 
case, C(H t ) = NIG (a, /3, St, fit) and the characteristic function resumes the 
form 



wnt. exp(5t v / a 2 - f3 2 ) 



<p Ht (u) = e mfMt ^ 7 v . (4.2) 

Y exp(5t^a 2 -(J3 + iu) 2 ) V ' 

Another interesting subclass is given by the hyperbolic distributions which 
arise for A = 1; the hyperbolic model has been introduced to finance by 



Eberlein and Keller (1995) 



Example 4.4 (CGMY model). Let H = (H t ) < t < T be a CGMY Levy 
process, cf. Carr, Geman, Madan, and Yor (2002) another name for this 
process is (generalized) tempered stable process (see e.g. ICont and Tankovl 
(|2004p ). The Levy measure of this process has the form 

X CGMY (dx) = C__l {a . >0} da: + C-_l {a . <0} dx, 

where the parameter space is C, G, M > and Y G (— oo, 2). Moreover, the 
characteristic function of Ht, t £ [0,T], is 

tpH t («) = exp (t C T(-y) [(M - iu) Y + {G + iu) Y - M Y - G Y ] ) , (4.3) 

for and the moment generating function exists for u G [-G, M]. 

The sample paths of the CGMY process have unbounded variation if 
Y G [1,2), bounded variation if Y G (0,1), and are of compound Poisson 
type if Y < 0. Moreover, the CGMY process is regular upwards if Y > 0; cf. 



Kyprianou and Loeffen (2005) Hence, by Statements 14.11 and 14. 2\ the laws 



of a CGMY Levy process, and its supremum, do not have atoms if Y G (0, 2). 
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The CGMY process contains the Variance Gamma process (cf. Madan 
and Seneta (|1990p ) as a subclass, for Y = 0. The characteristic function of 
H t ,t£ [0,T], is 



<p Ht (u) = exp [t C 



log 1 



lU 

M 



log 1 + 



lU 

G 



(4.4) 



and the moment generating function exists for u £ [—G,M]. The paths of 
the variance gamma process have bounded variation, infinite activity and 
are regular upwards. Thus, the laws of a VG Levy process and its supremum 
do not have atoms. 

Example 4.5 (Meixner model). Let H = (Ht)o<t<T be a Meixner pro- 
cess with C{H\) = Meixner (a, (3, 5), a > 0, — tt < (3 < n, 5 > 0, cf. 



tion of H t , t € [0,T], is 



Schoutens and Teugels (1998) and Schoutens (2002) The characteristic func- 



2St 



<PH t (u) 



COS ^ 



cosh : 



(4.5) 



and the moment generating function exists for u £ (— -, S+E\ The paths 
of a Meixner process have infinite variation. Hence the laws of a Meixner 
Levy process and its supremum do not have atoms. 



5. Applications in finance 

In this section, we derive valuation formulas for lookback options, one- 
touch options and equity default swaps, in models driven by Levy processes. 
We combine the results on the Wiener-Hopf factorization and the character- 
istic function of the supremum of a Levy process from this paper, with the 
results on Fourier transform valuation formulas derived in EGP. Note that 
the results presented in the sequel are valid for all the examples discussed 
in section UJ 

We model the price process of a financial asset S = (£t)o<t<T as an 
exponential Levy process, i.e. a stochastic process with representation 

S t = S e Lt , 0<t<T (5.1) 

(shortly: S = Soe L ). Every Levy process L, subject to Assumption (EM), 
has the canonical decomposition 

t 

L t = bt + ^~cW t + j J x(p-u)(ds,dx), (5.2) 

R 

where W = (Wt)o<t<T denotes a P-standard Brownian motion and \x de- 
notes the random measure associated with the jumps of L; cf. Jacod and 
Shiryaev IED031 Chapter II). 

Let M{P) denote the class of martingales on the stochastic basis ¥>. The 
martingale condition for an asset S is 

S = S e L G M{P) b + | + J (e x - 1 - x)X(dx) = 0; (5.3) 

R 
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cf. Eberlein et al. (2008) for the details. That is, throughout the rest of this 



paper, we will assume that P is a martingale measure for S. 

5.1. Lookback options. The results on the characteristic function of the 
supremum of a Levy process, cf. section [3l allow us to price lookback options 
in models driven by Levy processes using Fourier methods. Excluded are only 
compound Poisson processes. Assuming that the asset price evolves as an 
exponential Levy process, a fixed strike lookback call option with payoff 

(S T - K) + = (S e TT - K)+ (5.4) 

can be viewed as a call option where the driving process is the supremum 
of the underlying Levy processes L. Therefore, the price of a lookback call 
option is provided by the following result. 

Theorem 5.1. Let L be a Levy process that satisfies Assumption (EM). 
The price of a fixed strike lookback call option with payoff (|5.4p is given by 

Ct{S;K) = — / St iu ^ T (-u-iR)- — -— -du, (5.5) 

27r J T yiu — R)(l + iu — R) 



where 



( r 1 A [ e T(Y+iv) k(Y + zv,0) 

T A^oo 27T J Y + iv k(Y + iv,iu — R) 

-A 

for R £ (1, M) and Y > a*(M). 

Proof. We aim at applying Theorem 2.2 in EGP, hence we must check if 
conditions (C1)-(C3) (of EGP) are satisfied. Assumption (EM), coupled 
with Corollary 12.51 yields that (R) exists for R £ (— oo,M), hence 
condition (C2) is satisfied. Now, the Fourier transform of the payoff function 
f(x) = (e x - K)+ is 

f{u + iR) 



(iu - R)(l + iu — R)' 



and conditions (CI) and (C3) are satisfied for R € (l,oo); cf. Example 5.1 
in EGP. Further, the extended characteristic function ipj^ of Lt is provided 
by Theorem ES] and equals for R G (-oo, M) and Y > a*(M). Finally, 
Theorem 2.2 in EGP delivers the asserted valuation formula (15.51). □ 



Remark 5.2. Completely analogous formulas can be derived for the fixed 
strike lookback put option with payoff (K — S_ T ) + using the results for the 
infimum of a Levy process. Moreover, floating strike lookback options can 
be treated by the same formulas making use of the duality relationships 



proved in Eberlein and Papapantoleon (2005) and Eberlein, Papapantoleon, 



and Shiryaev (J2D08). 
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5.2. One-touch options. Analogously, we can derive valuation formulas 
for one-touch options in assets driven by Levy processes using Fourier trans- 
form methods; here, the exceptions are compound Poisson processes and 
non-regular upwards, finite variation, Levy processes. Assuming that the 
asset price evolves as an exponential Levy process, a one-touch call option 
with payoff 

can be valued as a digital call option where the driving process is the supre- 
mum of the underlying Levy process. 

Theorem 5.3. Let L be a Levy process with infinite variation, or a regular 
upwards process with infinite activity, that satisfies Assumption (EM). The 
price of a one-touch option with payoff (|5.7p is given by 

A 

i r tj-R-iu 

BC T (S;B) = lim — / S^ +lu ^-r (u-iR)— du (5.8) 

' v ' ' a^oo2tt J ^ LtK ' R + iu v ' 

-A 

= P(L T >log(B/S )), 
for R G (0, M) and Y > a*(M), where ^>j )T is given by (|5.6p . 

Proof. We will apply Theorem 2.7 in EGP, hence we must check conditions 
(D1)-(D2). As in the proof of Theorem 15.11 Assumption (EM) shows that 
condition (D2) is satisfied for R G (— oo,M), while Theorem 13.61 provides 
the characteristic function of Lt, given by (|5.6p . Example 5.2 in EGP yields 
that the Fourier transform of the payoff function f(x) = l{x>iogB} equals 

and condition (Dl) is satisfied for R G (0, oo). In addition, if the measure 
P-£ t is atomless, then the valuation function is continuous and has bounded 
variation. Now, by Statement 14. 2\ we know that the measure Pj^ is atom- 
less exactly when L has infinite variation, or has infinite activity and is 
regular upwards. Therefore, Theorem 2.7 in EGP applies, and results in the 
valuation formula (|5.8[) for the one-touch call option. □ 

Remark 5.4. Completely analogous valuation formulas can be derived for 
the digital put option with payoff 1{ St<B }. 

Remark 5.5. Summarizing the results of this paper and of EGP, when 
dealing with continuous payoff functions the valuation formulas can be ap- 
plied to all Levy processes. When dealing with discontinuous payoff func- 
tions, then the valuation formulas apply to most Levy processes apart from 
compound Poisson type processes without diffusion component, and finite 
variation Levy processes which are not regular upwards. This is true for both 
non-path-dependent as well as for path- dependent exotic options. 

Remark 5.6. Arguing analogously to Theorems 15.11 and 15.31 we can derive 
the price of options with a "general" payoff function f(Lr). For example, one 
could consider payoffs of the form [(St — K) + ] 2 or STl^g T>B y, cf. IRaiblel 
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>C= ^ N l ~ J?~" n , (5-10) 



(|2000^ Table 3.1) and Example 5.3 in EGP for the corresponding Fourier 
transforms. 

5.3. Equity default swaps. Equity default swaps were recently introduced 
in financial markets, and offer a link between equity and credit risk. The 
structure of an equity default swap imitates that of a credit default swap: 
the protection buyer pays a fixed premium in exchange for an insurance 
payment in case of 'default'. In this case 'default', also called the 'equity 
event', is defined as the first time the asset price process drops below a fixed 
barrier, typically 30% or 50% of the initial value So- 

Let us denote by tb the first passage time below the barrier level B, i.e. 

t b = inf{t > 0;S t < B}. 

The protection buyer pays a fixed premium denoted by K, at the dates 
Ti,T2, . . . , T/v = T, provided that default has not occurred, i.e. Tj < tb- In 
case of default, the protection seller makes the insurance payment C, which 
is typically 50% of the initial value. The premium K, is fixed such that the 
value of the equity default swap at inception is zero, hence we get 

CE[e~^l {TB < T} ] 

where r denotes the risk-free interest rate. 

Now, using that l{ Ts <t} = l{s t <B} which immediately translates into 

P(t b <t) = E[l {TB < t} ] = E[l { s t < B} ], (5.11) 

and that 

T 

E[e~"» 1 {TB < T} ] = J e- rt P TB (dt), 
o 

the quantities in (|5.10p can be calculated using the valuation formulas for 
one-touch options. 
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